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Theory of ac-Stark splitting in core-resonant Auger decay under strong x-ray fields 
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4 In this work we report the modification of the normal Auger line shape under the action of an 

5 intense x-ray radiation. Under strong Rabi-type coupling of the core, the Auger line profile develops 

6 into a doublet structure with an energy separation mainly determined by the relative strength of 

7 the Rabi coupling. In addition, we find that the charge resolved ion yields can be controlled by 

8 judicious choice of the x-ray frequency. 
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10 INTRODUCTION 

11 The interaction of an atomic system with a radiation 

12 field in the regime of x-rays will lead to its ionization. 

13 The most dominant process will be, first, the ejection 

14 of an inncr-shcU electron (photo-clcctron) with the ab- 

15 sorption of a photon followed either by an intra-atomic 

16 Auger and/or a fluorescence transition. For relatively 

17 light atomic systems, the dominant decay channel of the 

18 single-hole singly-charged system is through an (radia- 

19 tionless) Auger transition, designated as 'normal Auger' 

20 which is a manifestation of electron-electron interaction. 

21 One variation on this scenario is to promote an inner- 

22 shell electron to an excited bound state, often denoted 

23 as Resonant Auger State (RAS), which can decay either 

24 through an Auger transition or by the emission of an x- 

25 ray photon. This process was first reported by Brown 

26 [l[ and since then a large number of investigations have 

27 taken place (see for example and references there in) . 

28 Under excitation by the strong radiation fields, now 

29 available, from Free Electron Laser (FEL) sources [3- 

30 01, the situation becomes considerably different. Rele- 

31 vant studies in the context of strong laser fields have 

32 been reported quite early Q and in response to recent 

33 developments in the x-ray wavelength regime a number 

34 of theoretical and experimental works have already ap- 

35 peared 0, fiol - fisj . In the simplest situation, Rohringcr 

36 and Santra in Ref. [lo| have studied the single-photon 

37 excitation of a neon K-shell electron to a RAS by an x-ray 

38 field and a multipeak Auger Electron Spectrum (AES) is 

39 obtained for the fields they considered. In this work, 

41 we examine the AES and the ionic yields in the case 

42 where a normal Auger process takes place (as opposed 

43 to the RAS process) which involves photoionization of 

44 the neutral from the K-shcU, followed by an Auger decay 

45 of the singly charged hole-system to the doubly charged 

46 ion. We demonstrate the emergence of ac-Stark split- 

47 ting (also known as Autler-Townes splitting [l3|) of the 

48 Auger resonance, resulting from strong Rabi-coupling of 

49 the apposite states in the resulting Ne+^ ion. In con- 

50 trast to RAS line-shape modification [lol |. the unusual 

51 phenomenon here is that the ac-Stark splitting is mani- 

52 fested in the kinetic spectrum of the Auger electron while 




Figure 1. (Color online) Schematic figure of the dominant 
channels involved in the interaction of neon with an x-ray field 
of frequency ca. 908.06 eV, chosen to match the Ne"'"^(ls~^ — 
3p) transition energy. 



53 strong Rabi-coupling occurs to one of the K-shell elec- 

54 trons of the doubly-charged ion. This effect requires an 

55 explanation on the basis of a two-electron representation 

56 of the Auger-electron ejection, instead of an ambiguous 

57 one-electron picture implied by the Rabi-type oscillation 

58 of an inner-electron. A detailed study of the two-electron 

59 representation, in a different context and formalism but 

60 similar physical background, can be found in Ref. jl5|. 

61 We have chosen as the target of our study neutral neon 

62 in it's ground state \G) = \^ie{ls'^2s'^2p^ , ^So)) and con- 

63 sidcr radiation with a photon energy of approximately 

64 w = 908 eV (see Fig. [1}. The ionization potential for 

65 the creation of a K-shch hole |z) = |Ne+(ls2s22p^, ^P)) 

66 is i?^*^ = 870 eV. Note that all the energies of the neon 

67 states are given relative to the neutral neon ground state. 

68 Thus, the x-ray radiation will eject a K-shell photo- 

69 electron with a kinetic energy of ca. £i = 38 eV. lon- 

70 ization from the outer shells is also possible but with 

71 much lower probability The generated K-shell hole 

72 Ne^ will decay, by filling the Is vacancy, predominantly 

73 to the doubly charged neon state \a) = |Ne+^(ls^2s^2p^, 
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74 ^D2)) with energy E'^") = 65.35 eV, and Auger decay 

75 width r,a = 0.27 eV 0. In addition, the Ne+ ion with 
75 one K-shell hole can also decay through fluorescence by 

77 emitting x-ray radiation [l6j . This channel is about 55 

78 times weaker than the Auger decay rate. Let's designate 

79 the sum of all decay channels of the \i) state as Fj. Nor- 

80 mally, this Auger transition of the singly charged ion to 

81 the doubly charged ground state is accompanied by the 

82 ejection of an electron (Auger-elecron) with a kinetic en- 

83 ergy ei°' = i?'*) — ~ 804.65 eV and a Lorentzian line 

84 profile. In the present case we have chosen the photon 

85 energy to match the Ne~''^(ls~^ — 3p) transition energy. 

86 Around this energy there is a manifold of excited states 

87 which we denote collectively as \a'). We show in table 1 

88 those states which lie within a 1 eV band around the ex- 

89 citation photon energy. While the Auger state \i) decays 

90 with a lifetime of about 2.44 fs, the field, through a 

91 Rabi-type transition, creates a coherent superposition of 

92 the ionic ground state \a) and all accessible excited states 

93 \a'). This Rabi-oscillation of the K-shell electron between 

94 the bound states of the Ne"*"^ ion will induce an ac-Stark 

95 splitting manifested in the kinetic energy of the ejected 

96 Auger-electron. We should note that no post-collision in- 

97 teraction with the photo-electron is taken into account, 

98 as the energy of the photoelectron is too large for such 

99 an effect to contribute in the Auger's electron spectrum 

100 (and vice- versa). To complete the picture, the excited 

101 states I a') decay either through an additional Auger pro- 

102 cess (Fq/ ) or through further ionization by absorbing one 

103 more photon (7a')- Finally, direct creation of a double 

104 core-hole Ne+^ of the neutral, through photoabsorption, 

105 is not considered since the double K-shell ionization cn- 

106 ergy is 1863 eV [I^. 

107 We discuss below the manner in which the AES is mod- 

108 ified due to the Rabi-coupling of the Ne"'"^(ls~^ — 3p) 

109 states and study it's behaviour in a quantitative manner, 
no To facilitate the interpretation of the results, we note at 

111 this point that for a modification of the Auger line to ap- 

112 pear, many Rabi-oscillations should occur within the rel- 

113 evant Auger lifetime. Equivalently, in the energy domain, 

114 it is required that the energy separation of the Augcr- 

115 line splitting (roughly equal to Rabi coupling strength) 

116 should be larger than the Auger decay width or the x-ray 

117 bandwidth, whichever is larger. Of course the detailed 

118 properties of the system and the field will matter as well, 

119 however the rule of thumb, as expressed above will be in 

120 general true. 

121 The structure of the paper is as follows. In section II 

122 we present the theoretical formulation and develop the 

123 decription of the processes in terms of a time-dependent 

124 density matrix system of equations. In section III we 

125 show some of the results of our present study in the case 

126 of a single-mode coherent field. In the final section (Sec. 

127 IV) we summarize our findings and discuss, very briefiy, 

128 two aspects of the problem that are necessary to put the 

129 description of the problem closer to the actual expcri- 
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Table I. The above table lists the transitions from the \a) = 
Ne^^(^D2) ground state to its excited states Ne^^(ls~^ — 3p) 
around the photon frequency to = 908 eV. The fourth column 
shows the corresponding oscillator strengths. The data are 
calculated using the Cowan suites of codes [20| . 

130 mental conditions, namely the fluctuations present in a 

131 FEL field and the field's spatial dependance. 

132 



133 THEORETICAL FORMULATION OF THE 

134 DENSITY MATRIX EQUATIONS 

135 The density operator of the system is obtained in the 

136 basis of |G), |/), \A), \A!), \R), \Fi),i = 1,2 states. The 

137 state |G), with energy E'*-^^ represents the neon ground 

138 state. The state |/) = \i\ £,;), with energy Ei = iS^*) -I- e^, 

139 represents the K-shell hole Ne"*" (state \i) with energy 

140 i?(*)) and the photo-ejected electron \ei) having kinetic 

141 energy Ei. The state \A) = |a;e(i,eia), with energy 

142 Ea = E^"-^ -I- Ea + Eiai represents the Ne"*"^ ion in its 

143 ground state (state |a) with energy E'")) with an ejected 

144 Auger-electron of kinetic energy and the photoelec- 

145 tron having now kinetic energy Sia- Similarly, the state 

146 1^') = |a';ea',eia'), with energy Ea' = E'-"'^ + +Sia', 

147 represents the excited state of Ne"*"^ (state |a') with en- 

148 ergy E^'^ ^) with the Auger electron having kinetic en- 

149 ergy Ea' and the photoelectron having now kinetic en- 

150 ergy Sia'- It should be noted that in the definition of 

151 the photoelectron and Auger-electron states the appro- 

152 priate angular momentum quantum numbers, as they re- 

153 suit from electric dipole and Auger transition rules, are 

154 included. In addition, we also take into account the 

155 possibility of the involvement of further decay modes. 

156 In the present case, the K-shell hole Nc+ |z) can decay 

157 through fluorescence to Ne"'"(ls^2s^2p^), denoted here as 

158 \R). Moreover, the excited states \A') may further decay 

159 either through an Auger transition to Ne+'^(ls^2s^2p'^) 

160 denoted as \Fi) or through further photoionization to 

161 hollow K-sheU Nc+^{ls^2s'^2p^), denoted as jFa). The 

162 equations of motion for the density matrix elements are 

163 obtained from the Liouville equation ip(t) = [H(t),p(t)] 

164 with H{t) = H°+V+D{t), H° being the field-free Hamil- 

165 Ionian of neon, V the electron-electron interaction oper- 

166 ator and D{t) the x-ray field- atomic dipole interaction 
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operator. Inserting the above states into the Liouville 
equation we obtain: 

PGG{t) 2Ijn ^p GipiG, 

pii{t) = 2Im [DiGPGi] + 2ImYyiAPAi 

•A 

+ 2Im '^ DirPri 

PAA{t) = 2/m [VaiPia] + 2Im [Daa'PA'a] 
PA'A'it) = -2Im [Daa'Pa'a] + 2Im'^VA'FiPFiA' 

+ 2Im '^P A'F2PF2A' 

ipAA'it) = Eaa'Paa' + Daa'{pa'A' - Paa) + VaiPia' 

— ^^PAFiVfiA' - ^^PAFiDFiA' 

ipGi{t) = •■■ 



169 In the above expressions pkl , L = 

170 G, I, A, A' , R, Fi, F2 are the density matrix elements of 

171 the involved states while Dkl and Vkl represent electric 

172 dipole and Auger (intra-atomic) transitions between 

173 the states K,L, respectively. More specifically, the 

174 quantities Dgi,Via, Dir, Daa' ,Va'Fi, Da'F2 represent 

175 multielectron electric dipole (D) and Auger {V) transi- 

176 tion matrix elements. Within the present context we do 

177 not take into account any post-collision effects between 

178 the photo- and Auger-electrons as their contribution 

179 are expected to be negligible for the considered kinetic 
ISO energies. This assumption allows for a simplification 

181 of the transition matrix elements as for example for 

182 V/A = {i, ei\V\a, Ea, £ia) which reduces to Via = 

183 {i\V\a){eia\ei) = V^aS{eia - Along the same lines 

184 the dipole transition Daa' = {o-,£a,£ia\D\a',e'g^,£ia'), 

185 is approximated as Daa' = {0'\D\a'){ea\£a'){£ia\£ia') = 

186 ddd' 

5{£a — £a')5{£ia - £ia')- A detailed discussion of 

187 the dimensional reduction of these special kind of 

188 continuum-continuum matrix elements can be found in 

189 the appendix of Ref. |l5[- The summations involved 



here imply integration over the appropriate continua. 



As the total number of independent equations is 28, we 

192 do not present the explicit expressions for the evolution 

193 of the remaining density matrix elements as they are not 

194 essential at this stage. 

195 The density matrix equations, are a system of coupled 

196 integro-differential equations which are not amenable to 

197 an easy solution even by numerical means as it includes 

198 integration over multidimensional continua. It is thus 

199 our purpose here to transform the above system of equa- 

200 tions into a more tractable form. To this end, we adi- 

201 abatically eliminate the density matrix elements which 

202 are involved in the integrations over the respective con- 

203 tinua of the states. The procedure for adiabatically elim- 

204 inating these continua is a standard technique applied to 

205 describe the influence of a system with infinite degrees 

206 of freedom on to a system with a small number of de- 

207 gre es and appears in many different contexts (see e.g. 

208 [l8|). Here, the reduced system is the one described by 

209 |G), |/), 1^) and \A') while \Fi),i = 1, 2 represent the 

210 dissipative environment. Within the present context of 

211 atomic continua, some of the details can also be found in 

212 [1^ . To proceed further, the radiation field is expressed 

213 as E(i) = e(f(i)e*"* -f £*{t)e-"^'^)/2 with e its polariza- 

214 tion vector and we transform to slowly varying variables 

215 by defining aki = pKLe~^'^^*,n = 0,±1,±2,±3, where 

216 n is chosen so that nu has the closest possible value to 

217 Ek — El. With the latter transformation we remove 

218 from the coherences the fast oscillation part of their evo- 

219 lution due to the frequency of the field (this is justified 

220 since for a frequency of 1 keV the field period is of the 

221 order of 4 as and all other time scales set by photoion- 

222 ization and Auger widths are of the order of 1 fs '--^ 1000 

223 as). Given that the radiation is in the form of a pulse we 

224 have kept the slowly varying envelope £{t) which in ad- 

225 dition may describe the stochastic properties of the field 

226 under consideration. This is however a problem which re- 

227 quires special care and postpone its discussion for now. In 

228 the present case, we assume a fully coherent single-mode, 

229 Fourier transform-limited, field. Then, by employing the 

230 rotating wave approximation (RWA) and keeping only 

231 the terms proportional to the first-order of the electric 

232 field, after tedious but straightforward manipulation we 

233 end up to the following set of equations for the reduced 

234 density matrix elements 
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235 where Q,gi{ei^t) = {g\D\i, ei)£{t) and fiaa'(i) = daa'S(t). 259 been calculated before-hand. The approximations lead- 

236 With 7g(t) = 2tt J dei\^gi{ei,t)\'^ we denote the pho- 260 ing to this set of equations require careful examination 

237 toionization width of the neon ground states relative to 261 of the appropriate range of radiation intensities. One 

238 the Ne+ K-shell hole state while 7^' = Fa' + 7a' (0 ^''^ approximation is to assume that the relevant continua 

239 is the sum of the Auger decay width to states \Fi) (Fa') 263 are smooth around the energies of the dressed bound 

240 and the photoionization width to states IF2) (7a') of the 264 states within an energy range comparable with the Rabi- 

241 excited states \a'). In addition, Ti = Tia + "/r, is the sum 265 coupling matrix element. This requires that transitions 

242 of the decay width of the ionic state \i) through Auger de- 266 close to ionization thresholds should not be considered. 

243 cay to states \a) and through fluorescence to states \R). 26? In addition, in the derivation procedure, we have ignored 

244 The quantity Via represents the strength of the Auger 268 terms proportional to the second order of the field, such 

245 transition of the hole state \i) to the particular ionic 269 that fiaa' < 1- The latter approximation, given the ma- 

246 ground state \a). Therefore, the quantity is expressed 270 trix element daa' ^ 0.06 a.u., will restrict the range of 

247 as Via = 2i:\Via\'^ ■ Finally the energy differences in the 271 the intensities where the working equations are applica- 

248 above equations now include all the shifts associated with 272 ble below to 4 x 10^* W/cm^. 

249 the Auger and dipole couplings of the relevant states with 

contiminuv si&tes Eki = Ek + Sk- [El + Si), I ^g, I, a' ^he system of equations that have been derived 

25iwith5, = 5,a + 5„and5a'=5a'/, +5a'/2. Thc cxact - must, simultaneously, be numerically integrated for 

252 definition of the shifts and widths are as below: 275 all different photo-electron and Auger-electron ki- 

276 netic energies, so as to provide the populations for 

g _-29_^ f^^ I-PgjP 277 aii{ei),aaai£a, £i),o-a' a' {£a,£i) at infinite times. Since 



^ " 2 »"-^o J """^ Eo + uj — Ej + irj' 278 in our case we are only interested in the Auger-kinetic 

T-^ r l^/Ap ^™ energy spectrum and ionization yields regardless of the 

Sia ~ *~2~ ^ j|"Jo J '^^^ El — Ea + irj ' photoelectron, we must integrate the fi- 

^ . ID P populations over the photo-electron kinetic ener- 

Sir — i—^ = lim dEfi — — —, 282 gies and determine the following probabilities an = 



2 '7^0 J Ej - Er + irj' 

2 



^ 283 / deiaaa{£i),craa{£a) = / d£iaaa{£a, £i) , CTa' a' {£a) = 

g ^ ^ Tg' _ f \Va'Fi 284 J d£iaaia'{£ai£i)- An alternative and more economical 

2 »?^o J Ea' - Ep^ + ill 285 way of obtaining the same results is to derive a coarse- 



way of obtaining the same results is to derive a coarse- 
--y^, / \Da'F2\'^ ^"'^ grained version of the present equations for these re- 

Sa'f2 —'i— — 1™ J dEp^ + uj — Ep. + IT] ' ™ duced, averaged over the photoelectron energies, den- 

288 sity matrix elements. In addition, the reduced set of 

253 where use of the well-known formula lim,,^o 1/{x + irj) = 289 the density matrix equations is amenable to further ma- 

254 ■p(l/x) — iiT6{x) must be made to split up the integrals 290 nipulation as it allows for the derivation of analytical 

255 into their real and imaginary parts. 291 expressions for long pulses or their averaging for stochas- 

256 At this stage, a working set of equations are estab- 292 tic pulses. Thus, the new reduced set of equations, is ob- 

257 lishcd and their numerical solution is feasible, provided 293 tained by first setting all the derivatives of the coherences 

258 that all the dynamical parameters of the problem have 294 equal to zero [except CTaa'{£aT£i,'t)\ and then integrating 
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295 over the photoelectron energy e^. To demonstrate the 

296 reasoning of setting the derivatives of the coherences to 



297 zero, we work out the evolution equation of tTgiiei.t) co- 

298 herence [Eq. (fTf| ]. We integrate Eq. (fTf|) in an interval 

299 t,t + T with T « ^aa,,^g^,^i^ and wc obtain: 



,7s(T)+r, 



rt+T 



) / dt'<jg,{e„t') 



^gi{£i,T) 



rt+T 



dt'<Jgg{t'). 



300 The ionization width 7g(T), the dipole flgi{ei,T) and the 

301 Sg{T), Si{T) ac-Stark shifts that are included in the defi- 

302 nition of Egi (r) were removed from the integral as their 

303 value doesn't change much between t and t + r as a re- 

304 suit of the slowly varying trasformation of the variables. 

305 Since it will always be \agi{ei,t)\ << agg{t), we can ne- 



time average of <7g 



'gi 



{ei,t) as: 

Og,(£„T)/2 



Eg,{T)+Lo^i{^g{T)+T,)/2 



(r), (2) 



where agi{ei,T) = J^'^''dt'agi{ei,t')/T 



and crgg{T) 
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gleet the left hand side and obtain the 'coarse grained' 309 dt'agg{t')/T. Thus, by setting the derivative of the 
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310 



coherence to zero, effectively, leads to a coarse-grained 

311 value for the coherence which follows adiabatically the 

312 ground state population. At this stage, integrating Eq. 

313 ()lb|) over time (in an interval [t,t + r]) and the photo- 

314 electron energy we obtain: 
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Eg,{T)+UJ~l{jgiT)+T,)/2 



<Jgg{T) 



= -T^a^^{T) + '-fg{T)agg{T) , 



(3) 
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/ d£i(Tii{ei,T) and cr^ 



315 where Uiiir) 

316 j^^^ dt' aii{£i,t'). To evaluate the integral we have as- 

317 sumed that f2gi(ei,r) is smooth over an energy range 

318 equal to the radiation's bandwidth (far from resonance 

319 structures in the continuum or ionization thresholds). 

320 Then by expressing Eg, as Egiir) = E'^^i^ + S g{T) - E^''^ - 



321 
322 



Si{T) - Si we have Im J dei\ngi{ei, t)\'^ / {Egi{T) + uj - 
*(7s(r)+r,)/2) ^7r|173,(£, =£;W+5,;-£;(s)-5g,T)|2 = 

323 27g(T). 

324 As the derivation is quite long and detailed for the re- 

325 maining coherences, we give here only the final result for 

326 the reduced (coarse-grained) set of density matrix equa- 



327 tions. These reduced equations are obtained working 

328 along similar lines as for the derivation of Eq. ([3|). In 

329 this derivation, we ingore the photoionization of the ex- 

330 cited |a'), {ja' = 0) as for the photon energies and the in- 

331 tensities considered is expected to be much less than the 

332 Auger decay transition, represented by Ta' ■ Therefore, 

333 after setting all the derivatives of the coherences equal to 

334 zero [Eqns. JlJ, ([Ig]), ([ih|, Q], we solve for the co- 

335 herences and substitute their values into Eqns (fTc]). (jldp 

336 and (|lep . Then we integrate our equations over the pho- 

337 toelectron's kinetic energy and obtain the coarse-grained 

338 (also changing r — >■ t) set of density matrix equations. 
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d-iiit) 



-r,;(Ti. 



0-a'a'{£a,t) = -E^'CTa'a' + /m [fi*, fTaa' ] , 



i(7aa'{£a,t) = {Sa' - i^)(Ta 



(4a) 

(4b) 

(4c) 
(4d) 
(4e) 
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339 The dynamics of the process are governed by the ioniza- 340 tion width of the neutral target 7g(t), the core Rabi- 
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341 coupling fla'{t) the intra-atomic decay rates Fi^Ta', 

342 the Augcr-ficld induced couphngs i^^i{t) and the Auger 

343 Aa{t) = E, - (Ea+ea) and field Sait) = [E^+u) - E^- 

344 detunings. Note that for notational simplicity we de- 

345 note the core Rabi coupling as Q,a' = ^aa' = ^a'a- We 

346 have defined n^,{t) = 2|V^„| V(A^,i7a') with A±(t) = 
34. ta~et+ ilt{t)l2, et{t) = ei"' + T f^i';^(i)]/2 and 

348 7±(i) = r, + [Ta' ± n^j{t)]/2. The real quantities v!-^) 

349 and Cl^^} arc defined in terms of the generalized Rabi frc- 

350 qucncy: 

it) = + = y {5,, p. (5) 

351 The AES, at detection time, is obtained by adding 

352 the contributions from the ground \a) and excited 

353 \a') states of the doubly-ionized neon: S{ea) = 
'^j^^ ^, dt'&jj{ea,t'). The populations of the 

355 states \a) and \a') are obtained as Pjj{t) ~ 

356 dt J dei&jj{ea,t'), i~a,a'. Although we numeri- 

357 cally integrate the above system of density matrix equa- 



358 tions Eqns @ the Auger spectra can also be cast in an- 

359 alytical form to a very good approximation. For exam- 

360 pie, for a non-decaying excited states (Ta' = 0) the gcn- 

361 eralized Rabi frequency becomes a purely real quantity 

362 Vta' = \J S"^, + 4|i7a' p while the AES become independent 

363 on the coherence evolution aaa' ■ The analytical approxi- 

364 mation consists of considering a pulse of constant ampli- 

365 tude £(t) which turns the coarse-grained system to a sys- 

366 tem of ordinary differential equations with constant co- 

367 efficients. Then the Rabi frequency, ac-Stark level shifts 

368 and the ionization widths become independent on time 

369 and in the expression for the AES only the population of 

370 \i) is time-dependent: 

/ + 00 
dt[ t) + CTa'a'{£a,t)\ 

-oo 

p /' + 00 

= Im[{Aa + Sa' - - n-)] / dtauit). 

ill Solving the Eqns (|4a|) and (j4b| for aii{t) we find that 

372 J_^dtaii{t) = 1/Ti and after some algebra we end up to 

373 the following analytical expression for the AES: 



S{ea) = 



47r 



.(0) 



s„,+n„ 



4 



(6) 



374 We should note here, that the predictions of the analyt- 

375 ical expression differ from the numerical solution in that 

376 it doesn't include the transient effects of the physical pro- 

377 cess which are expected to occur at times of the order of 

378 l/Ej: 



379 RESULTS 

380 In Fig. 2 we show the effect of the field strength in 

381 the Auger spectra for oj = 908.06 eV and assume a pulse 

382 envelope £{t) = sin^(7rt/T), with T being the total pulse 

383 duration equal to about 20 times the Auger life time 

384 - 48.8 fs (FWHM = 24.4 fs). More specifically in ah 

385 the following calculations the field included 10754 cy- 

386 cles. The ionization width of the neutral is given by 

387 7g(t) = 4.375 X 10~^i?^(i) a.u.. In the calculation we 

388 only include the |a') = |Ne+^(ls~^ - 3p), ^F ) excited 

389 state [state 3 of table 1] and assume its decay width 

390 to be zero (Fq/ = 0). For this, on-resonance process 

391 (Sa' = 0), the Rabi interaction energy is found to be 

392 fia' = 0.061£(i) a.u.. According to Eq. for peak in- 

393 tensities of/o = l.Ox 10l^ l.Ox lO^^, 3.51 x lO^*^, l.Ox 10^^ 

394 and 3.51 x 10^ W/cm^ peak separations (of equal height) 

395 of about 0.28,0.88,1.66,2.79 and 5.25 eV, respectively, 

396 should be expected for the AES. Particularly, for the 



397 lowest intensity of /q = 1.0 x 10^^ the peak separation is 

398 comparable to the Auger decay width (0.27 eV), thus the 

399 separation is hardly seen. The latter values coincide with 

400 the full numerical solutions, shown in Fig 2. In general, 

401 according to the above formula, for on-resonance condi- 

402 tions and long pulses (compared to l/F^) a change of the 

403 splitting by a factor around 3.2 should be expected, for a 

404 change of the peak intensity by one order of magnitude. 

405 For off-resonance conditions we have two unequal peaks 

406 with energy separation again determined from the gencr- 

407 alized Rabi-frcquency Cla' and relative height determined 

408 from the field detuning A^' . 

409 Next we turn to the case where the decay chan- 

410 nels of the excited core states |a') are present. We 

411 have found that for the intensities and photon energy 

412 considered, further photoionization of \a') is a much 

413 weaker channel compared with a RAS transition to Nc+'^ 

414 [iO, im . We take the RAS width to be a large fraction 

415 of Ne+3(ls2p22p4) ^ Ne+3(ls22p22p3) decay width and 

416 assume Ta' = 0.156 eV [2l[. In this case the life time of 

417 these excited ionic states \a') is about 4.12 fs, compara- 

418 ble to the Auger decay life time of interest here (~ 2.45 

419 fs) but still much shorter than the pulse duration. Given 

420 that the Rabi couphng saturates the Ne+^(ls^^ — 3p) 

421 transition very quickly (in the sense that their popula- 

422 tions are almost equalized) it can be expected that a 
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Figure 2. (Color online) Variation in the AES as the peak 
intensity of the pulse increases. The field has total duration 
about 48.8 fs (FWHM = 24.4 fs) and the photon frequency 
chosen equal to 908.06 eV, while F^/ — 0. The inset label give 
the — . o 
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Figure 4. (Color online) Populations of the Ne"*"^ and 
Ne+2(ls-^ - 3p) states as a function of time peak intensity 
Id = 5.0 X 10^^ W/cm^. All other parameters as in Fig. 3. 
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Figure 3. (Color online). Ionization yields versus the peak 
intensity of the x-ray field. Field parameters as in Fig. 2 and 
r„/ = 0.156 eV. Solid curves refer to a; = 908.06 eV while the 
dashed ones to oj = 904.06 eV. 



423 large portion of the population will very quickly decay to 

424 the Ne+3 ion. This behaviour is shown in Fig. 3 where 

425 the final populations in Ne"*"^ and Ne+'^ are plotted as a 

426 function of the peak intensity of the applied x-ray field. 

427 Thus for on-resonance cases (solid lines) we see that the 

428 majority of the population goes into Ne"'"'^ ion for all in- 

429 tensities considered. On the other hand, if we choose to 

430 detune the FEL to a photon energy of 904.06 eV (dashed 

431 lines), efficient population of the excited ionic state \a') is 

432 prohibited. This also causes an effective increase in the 

433 magnitude of the generalized Rabi frequency according 

434 to Eq. ([5]). In that case the situation changes dramat- 

435 ically. The relative population ratio is reversed for low 

436 intensity fields (< 10^^ W/cm^), with the Ne+^ domi- 

437 nating up to higher intensities where the ratio starts to 

438 decline in favor of the Nc+^ yield. Thus, it appears that 

439 a careful combination of the intensity and the photon cn- 

440 ergy of the field can control the relative populations of 

441 the triply and doubly ionic species of neon. It is worth 



443 to note here that for the present long pulses and for rel- 

444 atively strong fields (> 3.5 x 10^^ W/cm^), the depletion 

445 of the neutral neon might be significant. For example, 

446 for Iq = 3.5 X 10^^ W/cm^ the remaining neutral neon 

447 is 0.718 while for ten times stronger field Iq = 3.5 x 10^^ 

448 W/cm^ decreases to 0.036. Furthermore, it should also 

449 be noted that ionic species higher than triply and doubly 

450 ionized neon are not expected to contribute significantly 

451 as two- or multi- photon absorption is an unlikely ion- 

452 ization channel, given their ionization potentials and the 

453 large photon energy. Finally, a small fraction of Ne+^ 

454 is produced through the fiuorescence of the metastable 

455 K-shell hole Ne'*"^ ion. 



In Fig. 4 we show the populations of the Ne"*"^ and 
Ne+^(ls~^ — 3p) excited states as a function of time for 
a field of peak intensity 1.0 x 10^^ W/cm^. We have 
considered two different frequencies as in the case of 
Fig. 3, which represent the on-resonance (908.06, eV- 
solid lines) and off-resonance (904.06, dashed-lines) con- 
ditions. From this figure, we can see that the popula- 
tions of Ne"*"^ and Ne"'"^(ls^^ — 3p) are quickly (almost) 
equalized, thus allowing for the efficient production of 
the Ne"'"'^ ion through the Auger decay of the excited 
Ne+2(ls- 
tron (Fa 

we see that the same populations evolve differently with 
the amount of Ne"*"^ that is produced being significantly 
smaller than the Nc+^ yield for a broad range of inten- 
sities. Of course, even in the detuned case, when the in- 
tensity becomes higher, the Rabi amplitude will increase 
accordingly and again the quick transfer from the ground 
Ne+2 to the excited Ne+2(ls"^ - 3p) state wih allow the 
efficient production of the Ne^'^ ion. In that case Ne"'"'^ 
yield will surpass Ne"*"^ yield. By inspection of Fig. 3, 
the intensities that this overtaking of Ne+'^ takes place 
are beyond 3.5 x 10+^'' W/cm^. 



^ — 3p) state with the ejection of one more elec- 
= 0.156 eV). With off-resonance conditions 
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CONCLUSION 

We have presented a theory of the Auger kinetic spec- 
tra and ionization yields based on the time-dependent 
density matrix theory which encapsulates all the essential 
dynamical parameters of the physical processes that are 
involved such as photo-ionization, photo-excitation and 
Auger transitions. We have examined the AES and the 
ionic yields in the case where an inner-shell photoioniza- 
tion takes place followed by an Auger decay of the singly 
charged hole-system. We have demonstrated the emer- 
gence of ac-Stark splitting of the Auger resonance, result- 
ing from strong Rabi- couplings. In addition, ionization 
yields have been calculated for a range of intensities. We 
show how to control the branching ratios of various ionic 
species by varying dynamical parameters of the system, 
such as Rabi coupling and detunings. The theory was 
applied to the case of K-shcll ionization of neon with the 
photon frequency chosen to match the energy differences 
between ionic Ne"*"^ ground and excited states. 

In the present study, we have put aside the issue of a 
field undergoing fluctuations suitable for the description 
of more realistic situations. A more realistic approach 
is to assume a field with an amplitude undergoing ran- 
dom fluctuations In general, the fluctuations of 
the field give rise to a nonzero bandwidth (beyond the 
Fourier bandwidth) and intensity fluctuations [22|, [23j . 
The main differences are that the field will excite a num- 
ber of (nearby) excited states (see for example the most 
important ones in the table 1), basically those that lie 
within the x-ray field bandwidth and that the width of 
the Auger spectrum will be effectively increased with the 
x-ray bandwidth. The development of the appropriate 
theory, capable to describe the field fluctuations, requires 
a detailed formulation which is beyond the purposes of 
the present study and it will be the subject of a future 
work. The essential outcome of the present work, namely 
the emergence of the ac-Stark splitting in the Auger spec- 
tra under the presence of strong ionic Rabi-couplings, will 
remain and the stochastic nature of the field will mainly 
affect its observability. 

Finally, keeping the focus on the essence of the issue of 
experimental observability, we refrain from showing the 
influence of the volume integration in the final results, 
but its precise contribution needs to be evaluated when 



it comes to the actual experimental conditions. 
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